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Operator product expansion on the lattice: analytic Wilson coefficients 



1. Introduction 

The computation of structure functions is one of the challenging problems in hadronic physics. 
It can be related to the operator product expansion (OPE) of the product of two conserved vector 
currents between states of particles as nucleons or mesons. The OPE allows for the factorization 
of two length scales: the short (high energy) scale is associated to c-number expansion coefficients 
(Wilson coefficients) whereas the long (low energy) scale is put into the matrix elements of local 
operators [jD]. In the case of deep inelastic scattering - the relevant physical process to extract 
structure functions - the short scale is defined by the large momentum transfer Q 2 = —q 2 , the long 
range scale is given by the hadronic mass P 2 = M\. 

In lattice simulations, structure functions are calculated via their moments ^ n {q 2 ). Their 
generic form as determined by the OPE can be written as 

Jt n {q 2 ) = CP (q/n) Ai 2) (/i) + C " 4)( f M) Ai 4) 00 + . . . , (1.1) 

T 

where A^ and A^ are reduced nucleon matrix elements of local operators of twist-two and four, 

(2) (4) 

respectively. C,, and C n are the corresponding Wilson coefficients, }X is the renormalisation scale. 
Usually, the Wilson coefficients are calculated perturbatively. On the other hand, the computation 

(k) 

of the hadronic matrix elements A n is intrinsically non-perturbative. In the last years lattice QCD 
has been proven to be a very promising tool for this purpose [0]. 

In order to relate the computed quantities to physical observables one has to evaluate the 
Wilson coefficients and the operator matrix elements in the same scheme. Furthermore, on the 
lattice operators can mix with other operators of the same quantum numbers but with lower mass 
dimensions. Both the coefficients of leading operators and the expectation values of higher twist 
operators are plagued by unphysical renormalon contributions [||]. In the complete OPE sum these 
contributions must cancel [Q] which is not the case if one combines perturbative Wilson coefficients 
and non-perturbative hadronic matrix elements. Therefore, if one is interested in contributions of 



higher twist operators to the moments of structure functions as given in (1.1) it is best to carry out 
the OPE completely non-perturbatively on the lattice. This results in the following form for the 
moments 

•Ml 



Jl n {q 2 ) = C)t>(a,q)A)?>{a) + ~" 7*' Ajf ; (a) + . .. , (1.2) 



q 1 

with a denoting the lattice spacing. However, A^\a) and C^\a,q) can be divergent due to mix- 
ing of the twist-four operator with its twist-two counterpart, i.e. A^\a) ~ 1/a 2 and C« (a,q) ~ 
1 / (aq) 2 . So both the Wilson coefficients and the matrix elements must be computed non-perturbatively. 

Starting point of the calculation is the OPE of the Compton scattering amplitude of a 
photon (fn(q)) off a nucleon target ( L/V(P))) with momentum P: 

%v = (N(P)\f^q)P v (q)\N(P)) (1.3) 
= E^v^Afo/M) WO \^)l,..,»„\N{P)). (1.4) 



The sum in ( |1.3[ ) runs over all possible operators ^ n of spin n. The index m distinguishes 



operators of the same spin. In order to have reasonable convergence properties of the OPE (so that 
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it makes sense to truncate the expansion at a particular operator dimension) we need small P, i.e. 
P 2 -C q 1 . Moreover, W^v should not suffer from large lattice artefacts, so both P and q must be 
small compared to I /a. These conditions result in the constraint 

P 2 < q 2 < {27l/a) 2 . (1.5) 



In order to determine the moments of structure functions both nucleon matrix elements and 
Wilson coefficients must be computed. It has turned out that it is essential to control lattice arte- 
facts in the non-pertubative evaluation of the Wilson coefficients. Therefore we have to satisfy the 
constraint ( |L~5l ) as strictly as possible, but in addition a subtraction of lattice artefacts based on 
a perturbative calculation of the Wilson coefficients has proven indispensable. In this paper we 
present first analytic results for perturbative Wilson coefficients with lattice artefacts included so 
that 0(a 2 ) corrections can be computed. Our work is part of a project initiated by the QCDSF col- 
laboration []|] to calculate the moments of structure functions on the lattice with the OPE method. 



2. Calculation 

Let us define the reduced Wilson coefficients c(q 2 ) as 

C(a,q)=c(q 2 )C B oRN(a,q), (2-1) 

where Cborn carries the index structure. On the lattice both the non-perturbative C(a,q) and 
CBORN{a,q) have corrections of 0{a 2 ). So writing 

CBORN(a,q) = Cborn + ( a< ?) 2 ^born + ■ • • (2.2) 
C(a,q) = c^\q 2 )cfl RN +(aq) 2 c^{q 2 )c'i ) ORN +... , (2.3) 

and taking the ratio gives 

(2) 

n ^BORN 

which yields the desired 0{a 2 ) correction term. 

Because the Wilson coefficients are independent of the target states we calculate the tree level 
Compton scattering amplitude W^ v (a,p,q) with off-shell quark states of momentum p. The calcu- 
lation is performed in momentum space using symbolic lattice Feynman rules. The local operators 



n„ i n tne expansion (1.4) are restricted to contain at most three covariant derivatives in agree- 
ment with the intended numerical computations. We use Wilson fermions with free parameter r 
leading to the following tower of local operators 

• Unpolarised case: yly, yy^ D v y, V D M D V W, ¥7^i D v D m D p y 

• Polarised case: yy^y 5 y, ya^ v D a y, y/y M y 5 D v D m y, yo^DaDpDx y 
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with Dji the lattice covariant left-right derivative 

D M =D M - D M (2.5) 

and a MV = l/2[y M ,y v ]. 

The amplitude W^ v (a,p,q) is expanded in powers of sm n (ap) (n < 3). Using the relation 
Dv~ sin(a/? v ) we can project out all possible operators and their respective coefficients - the corre- 
sponding Born Wilson coefficients CBORN{a,q)- In the final step we express all operators in terms 
of operators transforming irreducibly under the hypercubic group H(4) [0]. 

3. Results 

In this section we give the results for the OPE of ^ v associated to the amplitude W^ v 

W^{a,p,q) = (p\^v\p) (3.1) 

calculated between off-shell states of momentum p (for a review see e.g. [f7|]). Even when restricted 
to operators with at most three covariant derivatives, the general result for the expansion of =^ v is 
too lengthy to be presented in closed form - it has 880 terms. Except for the part without covariant 
derivatives we give the results only for a special choice of the momentum transfer q - the maximally 
symmetric case 

q=(f,f,f,f), s = wif, c = cos/. (3.2) 



If not shown explicitly we set a = 1 and r = 1 . The momentum transfer ( |3.2[ ) is chosen to present 
results in shortest form, other choices are possible just as well. Due to lack of space we will give in 
this paper the first two orders in the OPE only. The complete results will be published in a future 
paper [g]. All Wilson coefficients are determined up to an overall factor which depends on the 
normalisation of the corresponding operators. We adopt the normalisation as given in [Q]. 

3.1 No covariant derivative 

Here we present the lowest order of the operator product expansion in the most general 



form 



/o)/ w 9 o - 8arcos(a<7u/2) 2 „ 
^v(a,q)/e r = -arS^y/ly 8^ v \j/lY + 

v-, 2ar cos(a<7 u /2) 2 cos(afl T ) 

L TP 5 M v?lV+ 

8ar sin( aq u / 2) sin( aq v / 2) 

— — - \jfl\lf — 



Q 2 

2ar i cos(aq t ) sin(aqn/2) sin(aq v /2) 



V Q 2 

2ar cos(aqu/2) sin(aq v /2) sin(aq^ 



Q 2 

( no., Il\ sin (na„) 



2ar cos(aq v /2) sin(a^/2) sin(aq v ) 



a 

2ia cost aq u /2) cos( aq v /2)sin(aq (j , 

V, (3-3) 

T,CT 
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with 

e 2 = G 2 («^)=Lsin(^ T ) 2 + ^(L(l-cosK T ))) • (3-4) 

e 2 , denotes the quark-photon coupling. 

Let us give the results for the momentum choice (3.2). For the diagonal part (f) (choosing 
jU = 1 as an example) we find the decomposition 

^Cfl/4=- 6(3 ~l (1 ~ c W , e}=4, 2 +i6(i- c ) 2 . 0.5) 

At this level the Wilson coefficients can simply be read off (the notation for the irreducible 
representations is explained in [^]): 



operator 


representation 


Wilson coefficient a expansion 




T (1) 


6(3-c)(l-c) 3o/i 5 („f\2\ 
Qj. 2 I 1 12 WJ J 



Reinserting the a dependence, the expansion in powers of a is obtained by expanding the sine and 
cosine functions. The off-diagonal part also contains contributions for polarised structure functions 
(pi = 1 , V = 2 as an example): 

f 0) , w , 2(3-c)(l-c) i(l+c)s, 
3$' (/)/<$ = - — — '-yW+^-^iwy^y- WAVY) (3-6) 



The Wilson coefficients are obtained as 



operator 


representation 


Wilson coefficient 


a expansion 


yrly 

¥Y3Y5¥,-¥YaY5¥ 


T (1) 
T (4) 


2(3-c)(l-c) 
i{l+c)s 

e 2 f 


f(i -U a f) 2 ) 

ij{l-%{af) 2 ) 



3.2 One covariant derivative 

In this case the local operators = \j/Yn D v y and & T „ VW = Y a ^v D a Y contribute. For the 
diagonal part ^l\f) we find four independent functions bi(c,s) which build the corresponding 
Wilson coefficients 

b x = 4/(l-c)(74-126c + 63c 2 -9c 3 )/0 
b 2 = -4/(6-8c + 3c 2 > 2 /0 
b 3 =4i(4-3c)s 2 /Q} 

b 4 = _4/(l-c)(4-9c + 3c 2 )/0 (3.7) 

The following table presents the Wilson coefficients belonging to definite representations of oper- 
ators building (/). 
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operator 

3(^11 + ^22 + ^33 + ^44) 
3(^11 + ^22 - ^33 - ^44), -fe{0n - 022) 

^5(^12 + £21), 75(^13 + £31), 75(^14 + 041 ) 

^(^23 + 012), ^(^24 + £42), 73(^4 + ^43) 
-4(^12 " ^2l), -4(^13 " -U^14 " ^4l) 



repr. Wilson coeff. a expansion 



T (1) 


2(^1+3^3) 




I(«/) 2 ) 


T (3) 
T l 




\5ia 2 1 1 
16 v 1 


-l(af) 2 ) 


T (6) 




8V2^ 


l(af) 2 ) 


T (6) 
T 3 


V2b 3 


2V2/2^ 


-\(aff) 


T (6) 
T l 


h{b2~b,) 


-4>d 


-U*f) 2 



For the off-diagonal part 5j 2 (/) we have 12 combinations b{(c,s) 

b x = -4/(l-c)(33-52c + 24c 2 -3c 3 )/0 
fc 2 = 4/(7-9c + 3c 2 >70 

63 = -4/(3-2c)s 2 /0 

64 = 6/(l-c) 2 (l+c)(2-c)/£/ 
b 5 = 2/(l-c) 2 (4-9c + 3c 2 )/2} 
b 6 = 2/(l-c) 2 (l+c)(4-3c)/0 
6 7 = -4(l-c)(19-18c + 3c 2 )V2/ 
b & = 2(l-c)(14-15c + 3c 2 )V<2/ 
6 9 = \2(2 -c)s 3 /Q 4 f 

bio = ~4s 3 /Q 4 f 

b n = -2(l-c)(4-9c + 3c 2 )s/Q 4 

b n = -2(4-3c> 3 /2/ (3-8) 
The corresponding Wilson coefficients are shown in the last table. 



4. Summary 

In this paper we have presented first results for perturbative Wilson coefficients in lattice QCD 
needed to compute moments of structure functions. We have used Wilson fermions to calculate 
the off-shell Compton scattering amplitude. The expansion in the lattice spacing a up to 0(a 2 ) can 
be used to determine the lattice artefacts in the numerical determination of the Wilson coefficients. 
As can be seen from the tables there are Wilson coefficients which are at least of 0(a). They are 
associated to the operators xjfly and \jf°iiv Dm W which do not occur in the continuum OPE for 
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operator 



repr. 



Wilson coeff. 



a expansion 



p n + ^22 + ^33 + ^44) 



'33 ' 



^(^12 + ^21 



^(013 + 031), 75(^14 + ^41 
^(^23 + ^32), ^(^24 + 



71' 



'13 



V2 
? 34 + ^43 

l) 



'23 



14 " 



/ 2 4 - 



1{23} 



+ ■ 



1{34} 



+ 



? 2{13p 
3{14} 



1{24} + ^2(14} 



2{13}' 



'2{34} 



+ ■ 



f 3{24}) 



f{ 14}' 



_ 2{14} 
-^{24} 



73«22 " <33)'-7g«22 + <33 " 
75 Kl " ^33)' 72 «11 + ^33 " 2^44) 

Z Epe(i,2,3) sgn(p) < 23 , \ E pG (i,2,4) sgn 

7=«22 + <33+ " 



r 

124 



fll + ' 



233 " 



144^ 



^244 ) 



'3 
'3 
'3 
'3 
'1 

,(6) 
-1 

'2 

-2 

'2 

'2 

-1 

-1 

'i 4) 
'4 

,(4) 
'1 

,(4) 



bi+b 5 
bi-b 5 
V2b 2 

75(^3+^4) 
75(^3+^4) 
V2b 6 

j= 2 (bs-b 4 ) 

^(2b w + b 9 ) 

V6bi2 

—V6bg 

-V2b u 

V2(b u +b 7 ) 

V2b u 



2 (b 9 



7^(4*11 



5l()J 

■2fc 7 ) 



"4>(1 + 

-4^(1 + 
r(l- 



ll(^) 2 ) 

^/) 2 ) 
2n 



2V2/2 

47M 1+ ^ (a/)2) 

4^2/ ^ 
3y3a /1 
2V2/ U 

^(i-f^ 2 ) 



l(af) 2 ) 



Uaf) 2 ) 



2 t(1 



3/ 
2^3«(1 

- a 
2^3/ 



3W) 2 ) 



-!(«/) 
a- 3 



4 {aff) 
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